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ABSTRACT 

A semi-analytic model of cluster cooling flows is presented. The model assumes that episodic 
nuclear activity followed by radiative cooling without mass-dropout cycles the cluster gas 
between a relatively homogeneous, nearly isothermal post-outburst state and a cuspy config- 
uration in which a cooling catastrophe initiates the next nuclear outburst. Fitting the model to 
Chandra data for the Hydra cluster, a lower limit of 284 Myr until the next outburst of Hydra 
A is derived. Density, temperature and emission-measure profiles at several times prior to the 
cooling catastrophe are presented. It proves possible to fit the mass M(a) with entropy index 
Pp~7 less than a to a simple power-law form, which is almost invariant as the cluster cools. 
We show that radiative cooling automatically establishes this power-law form if the entropy 
index was constant throughout the cluster gas at some early epoch or after an AGN activity 
cycle. To high precision, the central value of a decreases linearly in time. The fraction of 
clusters in a magnitude-limited sample that have gas cooler than T is calculated, and is shown 
to be small for T = 2keV. Similarly, only 1 percent of clusters in such a sample contain gas 
with Pp~Y < 2keVcm 2 . Entropy production in shocks is shown to be small. The entropy that 
is radiated from the cluster can be replaced if a few percent of the cluster gas passes through 
bubbles heated during an outburst of the AGN. 
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1 INTRODUCTION 

Observations from the Chandra and XMM-Newton missions have 
shown that the intergalactic media of cooling-flow clusters is not 
a multiphase medium of the type required by the theory of cool- 
ing flows tha t has been widely accepted for over a decade (e.g. 
Fabiar 1994). Moreover, when these data are taken together with 



earlier radio and optical data, it is now clear that over the lifetimes 
of these systems very little gas has cooled to temperatures much 
below the virial temperature. Consequently, these systems cannot 
be the scenes of a quiescent steady inflow regulated by distributed 
'mass drop out' as that theory pos t ulated. In the light of the new ob- 
servations ( Bohringer et alj 200 1 ;|Molendi & Pizzolato||200l|; Mc - 
Namara et al. |200l[ ; [Mazzotta et alj 2002t [Vlatsushita et al| |2002| ), 
increasing numbers of workers (e.g. Bohringer et al. [2002| and ref- 
erences therein) are accepting the view that catastrophic cooling of 
intergalactic gas to very low temperatures takes place only at very 
small radii, and that on larger scales the intergalactic medium is pe- 
riodically reheated by an active galactic nucleus (AGN) that is fed 
by central mass dropout. This gene ral picture has been ar gued for 
several t imes over the last decade ( Tabor & Binned 19931; B inney 
& Tabor |1 995|; Einneyj |l99fj, |l999j |Ciotti & Ostrikerl |1997|; Bin 



ney |2002| ), but the physics involved is complex because the prob- 
lem is inherently three-dimensional and unsteady, and many details 



remain to be filled in now that clearer observations are becoming 
available. 

Recently, a number of authors have re-investigated the possi- 
bility that heat transported from the cluster outskirts to the centre by 
thermal conduct ion could contribute significantly to re-plenishing 
radiated energy (Ivbigt et al.][2002[ |Fabian et al] Z002|; Z akamska & 
Narayan |20(4 |R uszkowski & Begelman| |2002| ). The classical ar- 
gument against thermal conduction remains that there is a narrow 
range of conductivity in which conduction is numerically signif- 
icant bu^fifils^0 i £hmina^ethe_ra & 
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significant conductive heat input to dense regions now, are probably 
incompatible with the existence of cooling-flow clusters since they 
would have caused cluster gas to evaporate early in the universe 
(Loeb 2002). In this paper we model the effects of AGN alone, and 
neglect thermal conduction. 

Although the mechanism by which the AGN heats the IGM 
is unclear - possibilities include the impact on ambient thermal 



plasma of collimated outflows from the AGN (Heinz et al 



Kaiser & Alexander 



1999 



1998 



Binney 199S ) and inverse Compton scat- 



tering b y thermal plasma of hard photons from the AGN ( Ciotti & 



Ostriker ]l997| , [200 1| ) - features have been observed near the cen- 
tres of several clusters that are almost certainly bubbles of hot, 
low-density plasma that have been created by the AGN. Simula- 
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Fi gure 1. The functio n M(a) for the Hydra cluster estimated from the data 



of IDavid et al 



(2001). The units of o" are 10 31 cm 4 g 2 ' 3 s 2 and in these 
units Oq = 17. The dotted line shows equation (Ml) with e = 1.48. 



tions (Churazov et al 
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PuThs^TaT] 200 1 i ^ruggen*et*aT 2002 ) suggest that these rise 
through the cluster's gravitational potential well on a dynamical 
timescale. The details of how a rising bubble mixes with surround- 
ing gas and disperses its energy around the intracluster medium 
cannot be securely deduced from simulations, because they depend 
on what happens on very small scales at the edge of a bubble. The 
general idea that heating by AGN results in a kind of convection in 
the ICM is probably correct, however. The purpose of this note is to 
present a simple semi-analytical model of the life-cycle of a typical 
cluster that this process drives. 

Section 2 defines the model, fits its initial condition to Chan- 
dra data for the Hydra cluster, and shows its observable charac- 
teristics at a number of later epochs. Section 3 discusses the cre- 
ation of entropy during an outburst of the AGN and estimates the 
fraction of the cluster gas that an outburst processes through bub- 
bles. Section 4 calculates the distribution of cluster-centre temper- 
atures in a magnitude-limited sample of clusters and shows that 
temperatures below ~ 2keV will very rarely be detected. Section 5 
similarly calculates the distribution of the minimum values of the 
specific entropy that will be detected in a survey. Section 6 sums 
up. Throughout we assume a flat cosmology with £2^ = 0.7 and 
H = 65kms" 1 Mpc" 1 . 



2 THE MODEL 

Gas that is in stable hydrostatic equilibrium is stratified by specific 
entropy s, such that s never decreases as one moves up the gravi- 
tational potential. Heating by an AGN disturbs this equilibrium by 
producing bubbles of high-entropy gas near the bottom of the well. 
On a dynamical timescale each bubble rises through the well, shed- 
ding as it goes gas of the appropriate value of s to each isopoten- 
tial shell. Following this rearrangement, the gas is in approximate 
hydrostatic equilibrium, but slowly radiates entropy. After a finite 
time, the gas at the bottom of the well, which initially has the least 
entropy and the highest density, has radiated all its entropy. Its den- 
sity is then extremely large and a further outburst of the AGN is 
stimulated when some of this dense gas accretes onto the central 
massive black hole. This new outburst starts the next cycle of the 
system. 

If P denotes pressure, p gas density and y the ratio of the prin- 



cipal specific heats, then s = constant + ln(fj), where = Pp - "* 
will be called the 'entropy index' and k B is the Boltzmann constant. 
Let M(a) be the gas mass with entropy index less than a; Fig. [j] 
shows an estimate of M(d) for the Hydra cluster from the data of 
(David et al. 200 1 , hereafter D200 1 ) plotted such that a straight line 



corresponds to the functional form 
M(a) =A(a-a ) E . 



(1) 



Since the data fall near a line, we assume that, at least during some 
phases in the evolution of the cluster gas, M(o) can be adequately 
fitted by this functional form. For the moment this choice is mo- 
tivated by the observations of the Hydra cluster. We will show in 
Section 2.3 that, even if at some early epoch, e.g. after an AGN 
activity phase, a = const, throughout the cluster, a power-law dis- 
tribution of M(o) naturally arises due to the effects of radiative 
cooling. The fit involves three free parameters, A, £ and CTq. The 
lowest entropy density index found in the cluster gas is Oq at the 
cluster centre. It is therefore intimately related to the time that must 
elapse before the next outburst of the AGN. The total mass of the 
cluster gas is given by 



-A(a n 



-co) 



(2) 



where G ma x is the highest entropy index of the cluster gas at the 
largest radius. This relation defines the constant A. 

Given a functional form M(a) one can find the hydrostatic 
equilibrium configuration of the gas in a given spherical gravita- 
tional potential <i>(r) as follows. Eliminating p from the equation 
of hydrostatic equilibrium we have 



dP 
dr 



dr ' 



A second equation relating P and M is 



m a 2 , i(P 
= Anr p = A%r — 

dr F \o 



l/Y 



(3) 



(4) 



For a given ^(r) we obtain a model atmosphere by solving the 
coupled differential equations ^ and (Q) from r = with the ini- 
tial conditions M(0) = and P(0) = Pq, where Pq is a trial value. 
At some radius R out the gas mass reaches the value M tot . The cen- 
tral pressure Pq is adjusted until the pressure at R ont is equal to a 
specified value Poo. 

It is instructive to recast equation ^ in terms of the gas tem- 
perature, T. Assuming the cluster gas to be ideal, we can eliminate 
P from equation (^) and find 



dlnT dlnrj 



dr 



dr 



k B T dr ' 



(5) 



where /j is the mean molecular weight. Since the second term on the 
right side of this equation is intrinsically positive, the indications 
from observations that dT /dr > at most radii requires dlna/dr 
to be significantly greater than zero. Now 



din a 

dr 



dlno 9 

= 471 pr . 

dM K ' 



(6) 



so dlnT/dr can be positive near the centre only if dM/dlno van- 
ishes rapidly as M — > 0. For our initial assumption forM(o), equa- 
tion (pi), this requirement translates to 



M 
Mtot 



1,/e 



Y-l ( M \ / nm.y \ 



amax/Oo-1 \4%r 2 pJ \k B T J dr 



dt> 



(7) 



As r — * 0, d^/dr — > constant in a Hernquist model and M/r 2 p ~ r 
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Figure 2. The fit between the unevolved model (curve) and the density in 
the Hydra cluster inferred from X-ray observations. Data points and errors 
from D2001. 



Figure 3. Comparison of the initial model with temperature of the gas in 
the Hydra cluster inferr ed from X-r a y obs ervations. Data points show the 



temperatures derived by David et al 



(2001 1. Vertical error bars are the 90% 



confidence limits and horizontal bars indicate the region within the cluster 
over which the temperature was averaged. The solid line shows the model. 



for any plausible density profile. Hence the right side of this in- 
equality vanishes with r when T tends to a constant, and as r if 
T ~ r 5 near the centre. If the left side is to vanish more slowly, we 
require £ > 3/(1 — 8) where p ~ constant and thus M ~ r 3 . 



2.1 Comparison with Hydra cluster 

In this section we calculate model atmospheres for the gas in galaxy 
clusters. We adjust the model parameters to provide a fit to the 
density and temperature distributions of the Hydra cluster as in- 
ferred from X-ray observations (D2001). We do not attempt a for- 
mal fit of the data as this would require folding our model results 
with the telescope response. We are only interested in a rough con- 
straint on our model parameters to find reasonable values. Thus 
varying the model parameters and fitting by eye is sufficient. From 
Fig. |l] we find that for the Hydra cluster e ~ 1 .5 in equation (|l|) 
and (To = 17 x 10 31 cm 4 g~ 2 / 3 s~ 2 . Out to the limit of the data at 
Pout ~ 230kpc, the ratio o" max /o"o ~ 13.5. We assume that the grav- 
itational pot ential has the form of the potential of an NFW dark- 
matter halo (Navarro et al. 1996). Then there are four free parame- 



ters to determine by fitting the observational data: the gas pressure 
at the centre of the cluster, Pq, the gas pressure at R out , Poo, the scale 
length of the NFW profile, r s , and the overdensity of the dark matter 
halo compared to the critical density at the redshift of the cluster, 
5 C . 

Figs ^ and [3] show comparisons of the model with the data. 
The good fit is achieved for Py = 8.8 x 10~ 10 ergcm~ 3 and = 
2.6 x 10~ n ergcm~ 3 . The length scale of the NFW profile is r s = 
79.3 kpc and the overdensity of the dark matter halo 8 C = 8.5 x 
10 4 for our chosen cosmology. Both values are consistent with the 
findings of D2001. From these parameters we find the gas density 
in the cluster centre is po = 1.7x 10~ 25 gem -3 , which corresponds 
to an electron density of 0.085 cm -3 . The gas temperature at the 
centre is k^To = 3.4keV. The mass enclosed within P ollt = 230 kpc 
isMtot = 5.8 x lO 12 fc _1 M . 

Since in our model p tends to a constant at small r, the in- 
equality (Q) cannot be satisfied with £ = 1.5. The upturn in the tem- 
perature at small r is a consequence of this fact. 



2.2 Cooling 

If A(r) is the usual cooling function, the rate at which specific 
entropy is radiated is 



s(r,t) = ffc B - = jr 1 ", 

where the electron density is 

rc e (r)=0.85^. 

m p 



(8) 



(9) 



From Fig. 9-9 of Binney & Tremaine (1987) we calculate the cool- 
ing function A(T) appropriate to metallicity Z = O.4Z0 (D2001). 
Then by integrating equation (^) at each value of r for a small time 
interval Sf, we can update the function a(M) that is initially speci- 
fied by equation ([!]) 

A(r), 



o(M) 



3 k B T 



8f, 



(10) 



where the quantities are evaluated at the radius that at time / con- 
tains mass M, Assuming P«, = const., we then use equations (^) and 
(Q) to reconstruct the density profile at time t + 8?, and so on. Note 
that R out will decrease as the cluster gas radiates energy, and the 
material outside P ou t slightly compresses the cluster. Fig. ^ shows 
the resulting predictions for the evolution of the density and tem- 
perature profiles of the Hydra cluster. 

The value of the central temperature halves, and the central 
density doubles in roughly 240 Myr. The central temperature then 
drops rapidly to zero in the following 44 Myr. We assume that at 
this time the central AGN erupts and restores the density and tem- 
perature profiles to distributions similar to those currently mea- 
sured. 

Fig. ^| shows the time-evolution of the central entropy index 
O"0- To a very good approximation, it falls linearly in time between 
t = and t = 260 Myr. For temperatures above about 2.6 keV our 
cooling function is dominated by bremsstrahlung, i.e. A(T) oc yJT. 
In this regime equation M implies 



rj ' 



. p2/5 1/10 



(11) 



The very weak dependence of 60 on either Po or CT leads to the 
very nearly linear time dependence of Oq. Even for temperatures 
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Figure 4. Density and temperature profiles for the Hydra cluster at the 
present time (full curves), after lOOMyr (dot-dashed curves), 200Myr 
(short-dashed curves), and 284 Myr (long-dashed curves). 
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Figure 5. The time-evolution of the central entropy index o"q. 



below 2.6keV, A(J) and therefore do are not strong functions of 
temperature, so the approximately linear time-evolution of <5q holds 
at nearly all times. Of course, equation ([II]) also holds for material 
further out in the cluster. A linear relation between t and a can thus 
be used for an accurate estimate of the time to elapse before the 
next nuclear eruption. 

Figure ^ shows the cumulative mass of the cluster gas as a 
function of a — rjy, where <5o(t) is the current lower bound on the 
entropy index. Only at times later than about 230 Myr is there a sig- 
nificant deviation from the power-law form of equation (hj). In fact, 



0.01 r 



0.001 



0.0001 



10.0 



Figure 6. Cumulative mass of the cluster gas as a function of the entropy 
index. Entropy index is shown in units of the current minimum entropy 
index of the Hydra cluster. The long-dashed line shows the initial power- 
law of equation (|lj). The solid line shows M(a) at ( = 284 Myr. The short- 
dashed line represents the power-law giving the best fit to the solid line. 
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Figure 7. The emission measure of the cluster gas at four different times. 



even at t = 284 Myr a power-law with exponent £ = 1.8 approxi- 
mates M(< a) to within a factor 2. For gas with an entropy index 
greater than 2fJo a power-law provides an excellent fit throughout 
the cluster evolution. 

The emission measure distribution is defined to be 



E(t,T) = I dV8{T'-T)nl 



(12) 



where T is the temperature of the gas in the volume element dV. 
Figure ^ shows the run of E with T at four times. The evolution 
of E(t,T) is dominated by the growth of a tail towards low T due 
to radiative cooling. Gas at T < 1 keV appears only at the very end 
of the evolution. Even gas with T < 2keV appears only at times 
later than about 245 Myr. We will show in Section ^ that this result 
expl ains the lack of observa tions of very cold gas in galaxy clusters 
(e.g. Bohringer et al. 2002). The secondary peak at temperatures 



below the main peak is caused by the break in the density distribu- 
tion at around 30 kpc (see Figure 0). 



2.3 The entropy index-mass distribution 

Inspired by observational data for the Hydra cluster, we have 
adopted an initial entropy index-mass distribution, M(a), of the 
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Figure 8. Cumulative mass as a function of entropy index for the cluster 
with an initially constant entropy index throughout. Entropy index is shown 
in units of the current minimum entropy index of the Hydra cluster. The 
lines show the distribution at intervals of 1 Gyr from 1 Gyr (top) to lOGyr 
(bottom) after the start of the calculation. 



form given in equation ([lj). We have shown that during the sub- 
sequent cooling of the cluster gas this power-law form is preserved 
(see Fig. ^|). Thus, any cluster for which M(o) is given by equation 
([!]) at one point in its evolution, will retain this form of M(o) also 
at later times. It is natural to ask whether the relationship between 
M and a — Go evolves towards the power-law form if a cluster is 
initially configured in a different way. To answer this question rig- 
orously, one would have to follow a cluster's evolution from in- 
finitely many initial conditions. A configuration of special interest 
is the adiabatic one in which all the gas has the same entropy in- 
dex o"o : convection, whether it is caused by AGN activity or infall, 
drives the sas towards this configuration. 

Fig. |8| shows the results of using the iterative procedure de- 
scribed in the previous section to follow the radiative evolution 
of an initially adiabatic cluster atmosphere. The parameters are 
set to the same values as for our study of the Hydra cluster ex- 
cept that the initial entropy index of all the gas in the cluster is 
3.7 x 10 33 cm 4 g~ 2 / 3 s~ 2 , 22 times larger than the current minimum 
entropy index in Hydra, and that the initial pressure of the gas at the 
centre of the cluster is 1.7 x 10 -10 ergcm~ 3 , a fifth of the central 
pressure in Hydra. These values were chosen so that the cluster gas 
evolves into a state that roughly resembles the current state of the 
Hydra cluster. 

From top to bottom, the curves in Fig. ^ show the M(a) re- 
lation in the cluster at t = 1,2,..., 10 Gyr after the start of the cal- 
culation. All distributions are virtually indistinguishable from the 
power-law form of equation (|l|) with an appropriate choice of e. In 
fact, the power-law form of M(a) is effectively established as early 
as ~ 0.1 Gyr (not shown in the Figure). At 10 Gyr the minimum 
entropy index of the cluster gas is equal to the currently measured 
value of o"o in Hydra. The entropy index also spans a range sim- 
ilar to that of the gas in the Hydra cluster. Only the exponent of 
the power-law distribution, e = 3.4, is considerably steeper. Start- 
ing from an initially uniform entropy index distribution, the pre- 
cise value of £ that results at later times depends sensitively on the 
chosen initial conditions. Also, as before we make the assumption 
that the pressure at the outer boundary of the cluster, P«,, is con- 
stant during the entire evolution of the cluster. Considering the long 
timescales involved and the significant contraction of the outer ra- 
dius due to the cooling process over this time, it is unlikely that 
Poo really remains constant. Any evolution of the pressure at the 



outer boundary will influence the value of e at later times. A more 
detailed study of these effects is beyond the scope of this paper. 

Thus a cluster that initially contains gas with a uniform en- 
tropy index, almost instantaneously develops a power-law distribu- 
tion for M(a). Once the power-law is established, it persists, with 
changing exponent e, until the gas at the centre of the cluster cools 
to very low temperatures. These results suggests that the observed 
power-law distribution M(a) of the Hydra cluster would appear to 
be a generic rather than an exceptional configuration. 



3 THE CREATION AND LOSS OF ENTROPY 

We now estimate the rate at which the AGN creates entropy in the 
ICM and compare it with the rate at which entropy is radiated by 
the ICM. We idealize the blowing of a bubble as a process in which 
N\ particles of the ICM are irreversibly heated to a high tempera- 
ture 7j from the temperature T$ at the base of the cooling flow. The 
heated bubble expands as it is heated and does work on its environ- 
ment. The interface between the bubble and the ambient medium is 
Rayleigh-Taylor unstable so it will dissolve on a dynamical time. 
From the sharpness of the edges of observed bubbles it thus follows 
that they inflate supersonically. Since material just outsid e bubbles 
appears not to be on a high adiabat (Fabian et al. 2000), we infer 



that the expansion is only mildly supersonic. The entropy genera- 
tion by weak shocks is small (see appendix), and we neglect it. 

The entropy of a non-degenerate ideal monatomic gas of N 
particles that occupies volume V at temperature T is 



S = Nk B 



-In 



V (2mkfiT\ 3 / 2 
N \ W J 



(13) 



Hence, when a bubble containing N\ particles is heated at constant 
pressure from T3 to T\ , the entropy created is 

S m = ^ l k B ln(T l /T 3 ). (14) 

Additional entropy is created as the bubble rises and mixes 
in with the ambient medium. Since most of the cluster gas falls 
within quite a small range in temperature, we simplify the equa- 
tions by supposing that the bubble mixes with material that all has 
temperature T 2 . Let the bubble mix with material that contains N 2 
particles and assume for simplicity that the mixing occurs at con- 
stant pressure. Then from eq. (JT^) it is straightforward to show that 
the entropy of mixing is 



l+N 2 T 2 /(NiT l 
I+N2/N1 
l+NiTi/(N 2 T 2 ) 



+ ^n ^— — j (15) 

Adding eqs ( |l4| ) and ( |l5j ) we get the total entropy created by a 
bubble. As an illustrative example, if T\ = 1072, ?3 = Z2/IO and 
N\ = N2/IO, then 5 in fl = 1 1 .5k B N\ and 5 mix = 2. l4k B N\ . The rela- 
tive contribution to the total entropy production from S m i x increases 
with N 2 /N\\ for example with T\ = \QT 2 again but N 2 = lOOJVi, 
Si n fl is unchanged while S m i x = 3.2k B N[. 

We may estimate the fraction of the plasma that is processed 
through bubbles during an outburst by finding the time taken for gas 
of density N 2 /V 2 at temperature T 2 to radiate the entropy created in 
an outburst, where V 2 is the volume of the gas. The rate of entropy 
radiation is 

._ A(r 2 )(0.52iV 2 /V 2 ) 2 l/ 2 



T 2 



(16) 
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Figure 9. Time f coo i required to radiate entropy generated by processing a 
mass fraction N\ /N2 of the cluster plasma through bubbles heated to 10 
times the ambient temperature. 



where 0.52 is the conversion factor between total particle density 
N2/V2 and electron density. The time to radiate the entropy created 
by a bubble is thus 



where 



Ni 

N 2 

N 2 

H In 

Ni 













+ ln 




73. 





l+N 2 T 2 /(N 1 T l ) 



I+N2/N1 
NiT 1 /(N 2 T 2 y 



1+Nx/N 2 



(17) 



(18) 



0.55A(r 2 )(A' 2 /V 2 )- 

Fig. ^ shows t c as a function of bubble mass fraction N\ /N 2 for 
T\ = 10T 2 , T 3 = T 2 /W with N 2 /V 2 = 0.02cm" 3 and T 2 = 3 x 10 7 K 
as numbers typical of the central part of Hydra A (D2001). A bub- 
ble mass fraction ~ 1 percent yields cooling times of the expected 
order, 200 Myr. 



4 DISTRIBUTION OVER T M in 

In data taken with Chandra and XMM/Newton it is found that the 
plasma in any given cluster is confined t o temperatures T m ax > T > 
7min> where typically r m j n > 10 7 K (e.g. Bohringer et al. 



2002 



and 

references therein). In terms of the model discussed here, r max and 
T m [ n are the largest and smallest temperatures at which there is sig- 
nificant emission measure. We now assume that the evolutionary 
time t of clusters is uniformly distributed in the interval (0,f coo i), 
where f coo i is the time required for the central temperature of a 
freshly reheated cluster to cool to zero. With this assumption and 
the approximation that we can neglect the evolution in the comov- 
ing density of clusters, we can calculate the fraction of observed 
clusters for which T m [ n lies in the range (T + dT,T). We assume 
that the clusters are drawn at random from a flux-limited catalogue 
that was compiled with an instrument of limited spatial resolution. 
Hence, we focus exclusively on the total flux coming from plasma 
at a given temperature, and ignore issues related to the spatial dis- 
tribution of the emission. This simplification is probably fairly re- 
alistic for current data. 

First we determine the emission measure distribution of our 
model as shown in Fig. ^| Multiplying E(t,T) by the cooling func- 
tion A(r) and integrating over T we obtain the cluster's bolometric 
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Figure 10. The fraction of detected clusters in which plasma cooler than T 
can be detected for two flux limits: 5 m ; n = 5 x 10~'°ergcm~ 2 s~' (dotted 
curve) and 5 m j n = 5 x 10~ 13 ergcm~ 2 s _I (bottom full curve). The top full 
curve shows the fraction of clusters in the fainter sample in which plasma 
cooler than T would be detected in a pointed observation ten times more 
sensitive. 



luminosity L(f). Dividing L by 4%D(z) 2 , where D(z) is the l umi- 
nosity distance of an object at redshift z (Carroll et al. 1992] ), we 
obtain the flux S(z,t) of a cluster observed at evolutionary stage t : 



S(z,t) 



4jiD(z) 



ATA(T)E(t,T) 



(19) 



We are interested in the number of clusters with 5 > S m ; n , so we 
invert this relation to find the redshift z max (0 out to which the clus- 
ter has 5 > 5 m ; n at stage t. With our assumptions, the number of 
clusters detected at evolutionary stages in (t + dt,t) is proportional 
to V[zmax(0] df, where V(z) is the comoving volume to redshift z- 

Similarly, the flux from plasma cooler than T m [ n in a cluster at 
redshift z is 



S{T m [ n ,z,t) 



1 

4kD(z) 2 



dTA(T)E(t,T), 



(20) 



which on inversion yields the maximum redshift Zmax(7min> f ) ou t 
to which the cool plasma can be detected. Finally, the fraction on 



observed clusters in which plasma cooler than T m 
is 



/(rutin) = 



J^dtV[z m!lx (T min ,t 
J^dtV[z max (t)} 



can be detected 



(21) 



The lower two curves in Fig. [l(] show f{T m { n ) for two values of the 
limiting flux 5 m i n , namely 5 x 10~ 13 and 5 x 10~ I0 ergcm~ 2 s~' 
(dotted). For values of 5 m ; n smaller than about 10~" ergcm~ 2 s~' 
the distance to which even the warmest plasma in Hydra can be 
seen is not large cosmologically, and f{T m \ n ) becomes insensitive 
to 5 m i„. At lower flux limits, cosmological dimming becomes sig- 
nificant, especially for the warmest plasma, and the fraction of clus- 
ters in which cooler plasma can be detected rises slightly. 

The lower two curves in Fig. [u^ predict that only ~ 10~ 3 of 
clusters with bolometric luminosity in excess of S m j n will have a 
flux in excess of S m j n from gas cooler than T ~ 2keV. In prac- 
tice when we examine a catalogue of clusters discovered in a 
magnitude-limited survey, gas at T < 2keV will have been detected 
in a fraction of clusters that exceeds 10~ 3 because many clusters in 
the catalogue will have been the subject of more sensitive pointed 
observations subsequent to their discovery in the original survey. 
To quantify this effect, the upper two full curves in Fig. |h] show 
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Figure 11. As Fig. [uj but giving the fractions of clusters in which plasma 
with entropy index smaller than o" would be detectable. The units of a are as 
in Fig. |lj In these units a = 0.9S2kQTn e 2 ^ with T in keV and n e in cm~ 3 . 



the fraction of clusters in which gas at temperatures below T gives 
rise to a flux in excess of 0.5S m j n and 0.1S m ; n (top curve). From 
these curves we see that if the pointed observations are ten times 
as sensitive as the survey limit, the fraction of clusters showing gas 
at T < 2keV rises to 5.6 percent, while the fraction showing gas at 
T < 1 keV is still only 0.04 percent. 



5 DISTRIBUTION OVER rj 



Kaiser ( 1991 ) pointed out that there is a gross discrepancy between 
the predictions of gravity-driven hierarchical models of structure 
formation and the statistics of X-ray astronomy. In particular, the 
theory of hierarchical structure formation predicts that the comov- 
ing density of X-ray sources was much higher at redshifts z> 1 than 
is observed, and the cluster mass-X-ray luminosity relation is less 
steep than that observed. Kaiser pointed out that these conflicts 
would be eliminated if structure-formation theory were extended to 
include heating of the intergalactic medium (IGM) by massive stars 
and AGN prior to and at the epoch of galaxy formation. As Kaiser 
pointed out, one should think of such heating as ra ising the adiabat 
on which the IGM lies. j^loyd-Davies et aL (200C ) have attempted 
to determine the adiabat of the IGM by measuring the entropy index 
a in galaxy groups and clusters at a distance of 10 percent of each 

system's virial radius, well outside the region in which radiative 

2/3 

cooling is important. They define an entropy index s = k^T/n e 
with units of keVcm 2 . Since n/n e =21/11, the conversion to our 
definition of a is 1 keVcm 2 = 9.82 x 10 30 cm 4 g~ 2 / 3 s~ 2 . Their 'en- 
tropy floor' of i-o.l ~ lOOkeVcm 2 is equal to the entropy index of 
gas at a radius of 200 kpc in our model of the Hydra cluster. This 
distance from the cluster centre corresponds to 10 percent of the 
virial radius of Hydra. The agreement is not surprising considering 
that the Hydra cluster was included in the study of L loyd-Davies 
et al. ( |2000[ ). 

Lloyd-Davies et al. specifically excluded from consideration 
gas which they considered dense enough to have cooled signifi- 
cantly because they wished to probe the state of the IGM before 
clusters formed. In the context of our model of a cooling flow, it is 
interesting to ask whether they would have encountered an effec- 
tive entropy floor even if they had not excluded cooling gas and had 
simply measured Oo, the minimum entropy index at the cluster cen- 
tre. In Fig. [H] we show the fraction of clusters that have detectable 



emission from plasma with entropy index lower than a - the cal- 
culation of this fraction proceeds in close analogy with that of Sec- 
tion^. We see that in a pointed observation ten times more sensitive 
than the survey from which the targets were selected, only 1 per- 
cent of clusters would have detectable emission from plasma with 
a < 2 x 10 31 cm 4 g~ 2 / 3 s~ 2 ~ 2keVcm 2 . Thus an 'entropy floor' 
with regard to the minimum entropy index of gas in galaxy clusters 
should be detectable in X-ray observations of clusters that have suf- 
ficient spatial resol ution. This entropy floor i s considerably lower 
than that found by ^loyd -Davies et alj (EOOOl). It is determined by 
the very rapid cooling of the central regions of clusters, rather than 
any pre-heating of the gas before the collaps e of the cluster that 



may explain the entro py floor at larger radii ( Vbit & Bryan 200 1 
McCarthy et"aT]|2002| ). 



6 DISCUSSION 

In the Hydra cluster we find that the distribution of mass in entropy 
index a has a simple power-law form. We find that this form is 
to a good approximation preserved when the cluster gas cools in 
hydrostatic equilibrium. Furthermore, the activity of an AGN at the 
cluster centre drives convective flows in the cluster gas. Continued 
convetion would eventually lead to a cluster with gas of a constant 
entropy throughout its entire volume. We show that starting from 
such a constant entropy distribution, the power-law form observed 
in the Hydra cluster arises naturally as the cluster gas cools. Thus 
the discovery that M(a) for Hydra fits a power-law is no accident, 
but would have been discovered if we observed Hydra at an earlier 
or a later epoch. 

We find that the central value of the entropy index, Co, is a 
linear function of time until right up to the final cooling catastro- 
phe that provokes an outburst of the central galactic nucleus. These 
results enable one to calculate the dynamics of a cooling flow with 
remarkable ease. 

We estimate that in the absence of heat sources Hydra has 
about 280Myr to run before there is a central cooling catastrophe. 
This result places a lower bound on the time between nuclear out- 
bursts. The actual inter-outburst time is likely to be longer both 
because there is residual heating by supernovae/galactic winds and 
a low-level AGN, and because we have no reason to suppose that 
Hydra is in its immediate post-outburst state; it has probably been 
cooling for some time. 

In trying to understand a cooling flow it is useful to focus on 
the system's entropy at least as much as on its energy. Between 
outbursts the story as regards entropy is simple: it is being radiated 
at a rate that can be readily calculated from the X-ray brightness 
profile. The story regarding energy is much more complicated: the 
gas losses energy radiatively but recovers some through the work 
done by both the surrounding IGM and the intracluster gravitational 
field. Our premise is that the nuclear outburst that is provoked by 
each cooling catastrophe restores the cluster gas to essentially the 
same state that it had immediately after the previous outburst. This 
premise requires that the entropy created by irreversible processes 
during each outburst is equal to the entropy radiated between out- 
bursts. 

We assume that the restructuring of the cluster gas during an 
outburst is effected by bubbles of plasma. We calculate both the 
entropy created when a bubble is inflated, and that created as it 
mixes in with, and heats, the bulk of the intracluster gas. The rel- 
ative sizes of these two entropy sources is a weak function of the 
mass fraction that passes through bubbles during an outburst, but 
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is typically of order 3-5. If a few percent of the cluster gas passes 
through bubbles during an outburst, the time between outbursts is a 
few hundred Myr, as the observations suggest. 

When gas cools from temperatures S3keV at which 
bremsstrahlung is the dominant cooling process, it spends only a 
very small fraction of the total cooling time at T < IkeV. In the 
case of cluster gas this effect is magnified by the large density con- 
trast between the cluster centre and the half-mass radius. Conse- 
quently, not only does gas at T < 1 keV exist only for a very small 
fraction of the inter-outburst time, but it is confined to an extremely 
small fraction of the total volume. We have used our simple models 
of cooling flows to calculate the fraction of clusters in a magnitude- 
limited sample in which gas cooler than temperature T would be 
detectable. We find that without pointed observations gas cooler 
than 2keV would be found in only ~ 10~ 3 of clusters; with pointed 
observations ten times more sensitive than the discovery survey, 
such cool gas would be detected in ~ 4 percent of clusters. Thus, 
we should not be surprised to find that gas cooler than ~ 1 keV has 
not been detected in clusters observed by XMM and Chandra. 

Similarly, we predict that pointed observations of a com- 
plete sample of clusters would detect gas at entropies below o ~ 
2keVcm 2 in only ~ 1 percent of clusters because low-entropy gas 
appears only fleetingly and in small volumes. This lower limit on 
the entropy is at a value of a that is a factor ~ 50 lower than the 
'entropy floor' claimed for intergalactic gas from X-rays studies of 
low-mass clusters at radii well outside the region where radiative 
cooling is important. However, with the spectral and spatial res- 
olution of XMM and Chandra the lower entropy limit should be 
detectable in samples of galaxy clusters. 
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APPENDIX: ENTROPY GENERATION IN SHOCKS 

A major source of entropy production is shock fronts. Consider a 
perpendicular shock front in which an ideal monatomic gas with 
specific entropy si , specific internal energy u\ and speed vi is con- 
verted into material in which the same quantities are denoted 52. "2 
and \<2- Then conservation of mass, momentum and energy across 
the shock front yield 



si _ n / 5r 2 /3 + l/2 
.2 " r 2 ^5n/3 + l/2 



(22) 



where the r,- = Uf/vf are the two roots of the quadratic equation 

(2r/3 + l) 2 

„ k—rk = constant - (23) 
5r/3+ 1/2 

Fig. [l2| shows the changes in u and s across a shock front as a func- 
tion of the Mach number^] of the shock, M = (vi — V2) /c s . The very 
slow rise of the curve for S2A1 at small values of M implies that 
weak sho cks are very ineffe c tive e ntropy generators. This is also 
shown in Landau & Lifshitz (1987). The much steeper rise in the 



curve for 112/ U\ shows that weak shocks do thermalize energy ef- 
fectively, but since they do so without significant entropy increase, 
the original bulk kinetic energy can be largely recovered by subse- 
quent adiabatic expansion. 



Our definition of M differs significantly from a common one that involves 
the upstream speed relative to the shock. With our definition, in the limit of 
a sound wave M — > 0, rather than M — > 1 as in the case of the alternative 
definition. 



